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Clapping modes in unconventional superconductors. 

A.V. Balatsky a , P. Kumar b and J. R. Schrieffer c 
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We consider a superconducting state with a mixed symmetry order parameter components, e.g. d + is or d + idf 
with d' = d xy . We argue for the existence of the new orbital magnetization mode which corresponds to the 
oscillations of relative phase <f> between two components around an equilibrium value of (j> = f • It is similar to the 
so called "clapping" mode in superfluid 3 He — A. We estimate the frequency of this mode u>q(B, T) depending on 
the field and temperature for the specific case of magnetic field induced d! state. We find that this mode is tunable 
with an applied magnetic field with wo(-B, T) oc BAq, where Ao is the magnitude of the d-wave order parameter. 
We argue also that similar filed induced clapping mode should be present in an organic p-wave superconductors. 
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Clapping mode in High-Tc superconduc- 
tor. The order parameter in high-Tc supercon- 
ductors is widely believed jj| to be of a d x i_ y i 
symmetry. However more careful consideration 
indicates that the symmetry of the state in high- 
Tc might be lower in a number of cases. The 
symmetry allows the secondary components to be 
generated whenever there is a perturbing field. 
These secondary components (is or id xy ) gener- 
ation has been addressed in recent literature for 
the case of inhomogeneity due to wall scattering 
pi or due to vortex texture S . Similarly, the gen- 
erated id' component near magnetic impurity and 
spontaneous violation of time reversal symmetry 
with global d + id' has also been discussed [0. 
Another possibility for a global d + id' phase has 
been pointed out 0-1?]] where the external mag- 
netic field applied to two-dimensional d-wave su- 
perconductor generates id' . 

We present here J7J an excitation which consti- 
tutes a direct test of the induced component of 
the order parameter. This excitation is uniquely 
tied to the existence of the secondary out of phase 
component of the order parameter. Consider the 
most general situation of a complex order param- 
eter which can be generally written as Ao + iAi 
where Ao is the original d x 2_ y 2 component and 
Ai is the induced s or d xy component which is 



orthogonal to the initial d x i_ y i state. Define the 
global (Josephson) phase of the order parameter 
v and the relative phase cj> as: 



A = [|A | + ea*>(ty(r))|Ai|]earp(ti/(r)) 



(1) 



If the order parameter A is to be defined at the 
Fermi surface, then the functions Ai,i = 0,1 
implicitly contain the angular dependence. The 
global phase v{r) can be position dependent and 
even singular as is the case for the vortex configu- 
ration. We will focus on the relative phase 4>. Its 
dynamics by definition is related to the appear- 
ance of the secondary component Ai. The dis- 
persion of this mode has a gap similar to the Lar- 
mour frequency in case of spins. In the sense of a 
magnetic excitation, that can respond to a time 
dependent magnetic field in a resonant manner, 
this mode is also comparable to the longitudinal 
NMR in 3 He - A §] . The details of the dynamics 
are clearly model dependent. 

To be specific we will focus on the field induced 
Ai = id xy secondary component in the bulk. In 
this case quasiparticle spectrum is fully gapped 
in the field with the minimal gap vanishing at 
zero field. The d + id' state breakes time reversal 
symmetry and has a finite magnetic moment M z 
perpendicular to the superconducting plane H,f7|. 
The mode discussed above in this case will corre- 



spond to the longitudinal oscillations of the con- 
densate magnetic moment around its equilibrium 
value. Below we focus on a 2d superconductor 
at the fields H ~ H C 2{T) and therefore we will 
ignore small (0(H — H C 2)) difference between in- 
duction B and applied field H||z. 

We find the resonance frequency for the rela- 
tive phase oscillations of 0. It turns out to be a 
gapped propagating mode with: 



, M clapping mode 

I similar to 3He-A 



uj 2 (B,k) = u 2 (B,T) + s 2 j (B,T)k i k j , 

TjB 



MB) 



N(0) 



A (B,T) 



Sy = dijS 



(2) 



(3) 



Here s — (a + bB 2 )Af ) is the mode velocity, a, b 
are some constants and r\ is a constant, a mea- 
sure of the strength of the interaction, which we 
discuss in more detail below; N(0) is the Density 
of States (DOS) at the Fermi level; both i,j refer 
to in-plane coordinates x, y. This mode is a lon- 
gitudinal magnetization oscillation SM z (t,r) ex 
5M z exp(iu)(B, k)t—ikr) and is tunable by the ex- 
ternal field. There are, in effect, two consequences 
arising from a non-zero r\ which may be used to 
estimate its magnitude. In the presence of a sec- 
ondary order parameter Ai = id xy , the temper- 
ature dependent upper critical field -ff c2 (T) has 
an additional contribution which is quadratic in 
(1 — T/T c ). This results in an upward curvature 
with a scaling field which depends on rj. 

We now will prove the existence of the "clap- 
ping mode" for d + ia" state. The free energy has 
the standard form : 



F = Fn 



B 2 



(4) 



F = ^|A | 2 + ^|A | 4 + ^lAAo^-ASI (5) 



F^llA^ + ^IAA^A?! 



(6) 



(7) 



Fi n t — — [D x , D y }A()Al 

where An, Ai are d and d' components of the or- 
der parameter, an = ao(T/To — 1) , Tn is the 



;»(B) -l,BA (B.T) 



Figure 1. Schematic representation of the longi- 
tudinal magnetization oscillations caused by the 
oscillation of the relative phase angle between An 
and Ai with an equilibrium value ir/2. The fre- 
quency of this mode is linear in the B field and 
is also linear in the magnitude of the An. This 
mode is therefore tunable by the field and by tem- 
perature. We estimate tx>o{B,T) ~ A\{B,T) and 
is in the Kelvin range. 



ordering temperature for A . The corresponding 
dx > = No is always positive, as are &n > and 
K l3 ,K' %J > 0. A = di - i^-Ai, with B = V x A. 
The interaction term in this form has been pro- 
posed earlier [|J[i]]. The interaction term, using 
[D x , D y ] = ieB can be written as 

F int = teB^AoAl + h.c. (8) 

which corresponds to a coupling of the magnetic 
field with an intrinsic orbital moment along z- 
axis:< M z >= ^A A^ + h.c.. 

Here we will assume that the amplitude for the 
secondary Ai has been developed and we will look 
at the relative phase <f> oscillations only. It is con- 
venient to introduce the respective phases of each 
component: 

An = \A \exp(i<f> ), A x = \Ai\exp(i<f>i) 

which are related to the introduced above v = 
4>o,4> — 0i — 0o- Derivation of the Eq.pft then 
proceeds as follows. We use the Josephson rela- 
tions for each component: 



dF 



V-i-, N t 



dF 



-,1 = 0,1 



(9) 



Where \ii are chemical potentials for particles in 
An,Ai condensate and similarly, N t are conju- 
gated number of particles. In the equilibrium, 



Fj nt = T) Aj.A. B sin<|) 



Figure 2. The fj n f profile as a function of the rel- 
ative phase angle <fi between Ao and Ai is shown. 
Although the minimum is reached at <fi = tt the 
hnite stiffness for <j> oscillations leads to the finite 
frequency mode at loq(B,T). We have ignored 
the dependence on the sign of B in the figure. 



where /io = /ii = fi, wc find for a relative phase 
motion: 



d/xo djM_ _ d/i dF 
[ dN 8^ dN ± l d<f> 



(10) 



where in the above we have taken into account 
the fact that although /io = \i\ their derivatives 
are different. The term in the brackets is on the 
order of N(0): [f^ + f^-2^] = p-i ~ JV(0). 
For the general values of the relative phase </> the 
terms in free energy F that contribute to the 
derivative in Eq.(|l0|) are Fi nt , Eq.(Q) and terms 
with gradients Kij , K'^ . 



F int = vB z \A \\Ai\sm l 



(11) 



with minimum F reached at (f> = —ir/2sgn(B z ). 
In general, there is an additional term propor- 
tional to the squares of the two (the bare and 
the induced) order parameters. However it con- 
tributes little new to the physics of the problem. 
Its effects are largely quantitative. We stress here 
that the phase </>o of Ao is not constant in the 
external field in the mixed state. Nevertheless 
the Fi n t m Eq. Jll] ) depends on the relative phase 
only. 

We find: 



= -? ? |Ao||A 1 |B z cos0-s 2 V 2 (; 
P 



(12) 



Here s is given by Eq.(0). Minimizing the free 
energy Eq.(QJ^J^) with respect to the magnitude 



from this equation, Eq.(g) follows immediately. 
The velocity s is field dependent and also can be 
used in experiments to detect the presence of id! 
component. 

We can estimate the magnitude of the energy 
gap up to a O(l) prefactors: 



u {B,T)~\A x {B,T)\ 



r/B 



A (B,T) (13) 



as one can easily see from minimization of the 
total F with respect to Ai. For estimate for the 
coupling constant r\ and magnitude of Ai see M. 
We estimate thus 



A (B,T) - [ e )( H c2 } 



(14) 



This result turns out to be very similar to the 
"clapping mode" in 3 He— A, where frequency was 
found to scale with the gap in the whole temper- 
ature range as well || . We also find asymptotics 
in H: 



LOq 



H/H c2 

\/H c 2 



_ H<H c2 

H H — ► H C 2 



(15) 



and we note immediately that the mode frequency 
ujq <C Ao Therefore this mode will be sharp. The 
damping coming from the low-lying quasiparticles 
in the nodes of d-wave gap will not affect this 
mode because the phase space for the decay will 
be small. 

We assumed that once the field induced com- 
ponent is present in GL region it will persist to a 
lower temperatures. We present Eq.(n5|) as an or- 
der of magnitude estimate only. From Eq. dlj) it 
follows that the gap a->o can be made in the range 
of 0.1 - 1K(2 - 20GHz) in the field H = 1 - 10T 
similar to the Ai estimates fft. 

The longitudinal oscillations of magnetization 
M z will also lead to the resonance at u>o(B,T) 
in the AC susceptibility of the superconduct- 
ing state. Experimentally the proposed clapping 
mode can be observed in NMR in the field. For 
any of the resonance techniques used to search 
for the "clapping mode" resonance at wq(B,T) 
Eq.(|lJ) the vortex cores contribution would be 
the biggest source of background. 

Clapping mode in p-wave superconduc- 
tor. Here we conder the case of quasi 2- 



dimensional p-wave superconductor in an exter- 
nal field. We identify here Ao oc p x and Ai oc p y . 
Assume that the zero field state has real order 
parameter that transforms as Ao oc p x . We 
argue then that external magnetic field will in- 
duce secondary component Ai = ip y |l(J. The 
state Ao + Ai will have a finite magnetic mo- 
ment < M z > that can couple to magnetic field. 
The free energy term driving the secondary com- 
ponent is given by Eq.(g) and we find Ai ex -BA 
in case of p-wave superconductor. The results 
for the clapping mode are similar to the case of 
d-wave superconductor and will be presented in 
more details elsewhere (HI. We note that the 
similar clapping mode for the p-wave state that 
violates time reversal in zero field, e.g. p x + ip y , 
was considered in [ pd| . 

In conclusion, we considered the relative phase 
oscillation mode that can be used to detect the 
secondary component in the time reversal violat- 
ing superconducting state, such as d + id' or 
p + ip' . Two components of the order parameter 
are characterized by respective amplitudes and 
phases. The relative phase between two compo- 
nents can oscillate around its equilibrium value 
<j> = ±-|. This mode is very similar to the "clap- 
ping mode" in superfluid 3 He — A. For a specific 
model we choose the case of field induced id' (for 
singlet pairing) or ip' (for triplet) component. We 
show how the relative phase mode frequency is 
governed by external magnetic field and d-wave 
gap magnitude luq(B,T) oc -BAo and is there- 
fore tunable by external magnetic field. Among 
other probes this mode could be experimentally 
detected by investigating the ac magnetic suscep- 
tibility and possibly by ultrasound attenuation in 
the in the mixed state as a function of applied 
field H. 

Work at Los Alamos was sponsored by the US 
DOE. 
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